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Mantle  Shear  Wave  Velocities  Determined  from  Oceanic 
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Palisades,  New  York 

Yasuo  Sato 
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Abstract.  Dispersion  of  Love  and  Rayleigh  waves  liZifeiBu  computed  for  a  model  of  the 
suboceanic  crust  and  mantle.  An  important  feature  of  this  model  is  a  low-velocity  layer  in  the 
upper  mantle.  These  calculations,  which  include  the  effect  of  ephencity,  are  consi.4ent  with 
oceanic  Love  wave  dispersion  data  for  periods  from  10  and  800  seconds  to  withia  a  few 
hundredths  of  a  kilometer  per  second.  The  effect  of  sphericity  on  oceanic  Love  waves  is  of  great 
importance  for  periods  as  short  as  10  seconds  because  of  the  penetration  of  these  waves  into  the 
low-velocity  channel.  Computations  for  the  same  model,  which  consists  of  homogeneous  iso¬ 
tropic  layers,  are  also  in  accord  with  the  oceanic  Rayleigh  wave  data  for  periods  from  30 
to  150  seconds  to  within  a  few  hundredths  of  a  kilometer  per  secondj^he  results  can  be 
further  improved  by  minor  modifications  of  the  model.  Since  the  calculaftons  for  a  spherical 
earth  are  consistent  with  oceanic  Love  and  Rayleigh  wave  data,  apparent  l^crepancics  which 
result  from  computations  for  flat  layers  are  resolved  without  recourse  to  V  difference  in  SH 
and  SV  velocities  in  the  suboceanic  mantle.  R.ecent  calculations  and  observ^ons  of  Love  and 
Rayleigh  wave  dispersion  by  Brune  and  Dorman  for  the  mantle  beneath  theXCanadian  shield 
.■3:;..werl  (1)  that  a  low-velocity  zone  is  required  and  (2)  that  no  difference  bet^en  SII  and  SV 
velocities  is  neeaed.  Thus  their  study  and  this  paper  show  that  these  two  conclusions  are  appli¬ 
cable  to  large  areas  of  the  world.  The  wave  guide  of  the  G  wave  for  periods  from  60  to  300 
seconds  is  the  free  surface  and  the  steep  gradient  in  .shear  velocity  in  the  upper  half  of  the 
mantle.  At  shorter  periods  the  fundamental  oceanic  Love  mode  exhibits  large  particle  ampli¬ 
tudes  in  the  low-velocity  channel.  A  comparison  of  the  upper  mantle  beneath  oceans  and  con¬ 
tinents  confirms  previous  results  which  indicate  that  channel  velocities  are  smaller  beneath  the 
oceans. 


Introduction.  The  physical  properties  of  the 
cnist  and  rnantic  beneath  oceanic  areas  can  be 
iiive.stigated  with,  the  aid  of  Love  .and  Rayleigh 
waves.  In  this  study  theoretical  phase  and 
group  velocities  for  a  model  v/hich  contains  a 
low-vc-Iocity  layer  are  compared  with  experi¬ 
mental  data  for  a  broad  range  of  periods.  The 
effect  of  sphericity  on  Love  and  Rayleigii  waves 
has  been  taken  into  account  in  the  theoretical 
calculations. 

In  this  paper  considerable  emphasis  has 
been  placed  on  the  Love  wave  solution  for 
a  spherical  earth  as  compared  with  the  result 
for  fiat  layers.  The  effect  of  sphericity  on 
the  phase  and  group  velocities  of  oceanic  Love 
waves  is  not  negligible,  even  for  periods  as  short 
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as  10  seconds,  because  of  the  penetrjition  of  the 
wave  disturbance  into  the  low-velocity  channel. 
In  fact,  sphericity  can  resolve  the  apparent  dis¬ 
crepancies  discussed  by  Ewing  and  Landisman 
[1961]  which  result  from  computations  of  Love 
and  Rayleigh  wave  dispersion  for  flat  layers. 
When  sphericity  i.s  taken  into  account  there  is 
no  longer  a  necessity  for  a  difference  in  SH  and 
velocities  in  the  suboceanic  mantle. 

Recently  Brune  and  Dorman  [1962],  using 
high-precision  Rayleigh  and  Love  wave  data, 
found  that  no  difference  between  SH  and  SF 
velocities  is  needed  for  the  mantle  beneath  the 
Canadian  shield.  They  also  demonstrated  that 
a  low-velocity  layer  is  required  for  this  region. 
Thus  their  study  and  this  paper  show  that  these 
two  conclusions  are  applicable  to  large  areas  of 
the  world. 

For  the  caleulation.s  ii.scd  in  this  study,  the 


5257 


r)2o8 


SYKES.  I.ANDISMAN,  AND  SATO 


ilistribiition  of  vclocitiP8  and  densities  for  the 
suboceanic  cnist  and  mantle  was  selected  with 
the  aid  of  several  different  criteria.  The  physical 
parameters  for  the  sr  Mmentary  and  crustal  lay¬ 
ers  of  case  122  were  coosen  to  be  in  accord  wdth 
the  results  of  seismic  refraction  work  at  sea 
and  with  the  short-pt'riod  oceanic  stirface  wave 
results  of  Oliver  and  Dorman  [1961].  For  the 
mantle  the  distribution  of  physical  parameters 
for  case  122,  which  is  consistent  with  data  for 
teleseismic  body  waves,  is  the  same  a.s  that  for 
oceanic  ca.‘;c  8099  of  Dorman  et  al.  [  1960] .  These 
authors  showetl  that  case  8099  was  in  accord 
with  the  group  velocity  flata  for  mantle  Ray¬ 
leigh  waves  determined  by  Ewing  and  Press 
r  1954a,  6]  and  with  the  oceanic  Rayleigh  wave 
observations  of  Sutton  et  al.  [1959].  One  of  the 
most  important  result^!  demonstrated  by  Dor¬ 
man  et  al.  was  the  requirement  of  a  low-velocity 
channel  for  the  mantle  beneath  continents  and 
oceans. 

Recent  obser\'ational  work  has  produced  ad¬ 
ditional  Rayleigh  wave  data  for  the  oceans,  some 
of  which  includes  phase  velocities  avS  well  as 
gr  up  velocities  [Aki,  1960a,  b;  Aki  and  Press, 
1901,  Kuo  ei  al.,  1952].  Ki:o  ct  al.  have  shown 
that  case  8099  explains  the  major  features  of 
oceanic  mantle  Rayleigh  wave  dispersion.  The 
present  study  indicates  that  theoretical  Love 
and  R.aylcigh  wave  phase  and  group  velocities 
computed  for  case  122  sati.sfy  available  oceanic 
dispersion  data.. 

Long-period  Love  waves,  especially  G  waves, 
have  been  studied  since  1920  by  Gutenberg,  By- 
erly,  Imamura,  Wilson,  Sato,  and  others.  Love 
wave  calculations  by  Sato  [1958],  Landisman 
and  Sato  [1958],  and  Landisman  ct  al.  [1959] 
I'.ave  contributed  to  the  present  knowledge  of 
the  upper  ma,ntlc.  With  tlu'  exception  of  tliose 
of  Sato  [1958],  these  Love  wave  calc.ilations, 
studies  of  teleseismic  bofly  wave  times  and  am¬ 
plitudes,  and  the  Rayleigh  wave  computations 
cited  above  have  sliown  that  beneath  the  crjst 
the  uiqKT  mantle  has  high-velocity  material 
which  overlies  a  region  of  lower  velocities.  At 
even  greater  depths  higher  velocities  are  also 
required.  This  distribution  of  velocities  is  com¬ 
patible  with  those  inferred  from  body  wave 
studies  by  Gutenberg  [1953]  and  Lehmann 
[1955],  A  general  conclusion  to  be  drawn  from 
these  .surface  wave  investigations  is  that  the 
channel  velocities  in  the  upper  mantle  are 


•smaller  beneath  the  oceans  than  beneath  the 
continents. 

The  first  studies  in  which  observational  group 
velocity  data  were  compared  with  calculations 
for  flat  layers  were  restricted  either  to  Love 
waves  or  to  Rayleigh  waves.  Most  of  the  result¬ 
ing  independent  velocity  distributions  con¬ 
firmed  the  existence  of  the  low-velocity  layer, 
but  they  were  quite  dissimilar  in  detail.  For  ex¬ 
ample,  Dorman  et  al.  [1960]  demonstrated  that 
case  38  km-XII,  which  had  been  showm  by 
Landisman  and  Sato  [1958]  to  explain  conti¬ 
nental  Love  wave  data,  was  not  in  accord  with 
continental  Rayleigh  wave  data.  Ewing  and 
Landisman  [1961]  discussed  this  problem  fur¬ 
ther,  paying  particular  attention  to  the  mantle 
structure  beneath  oceans,  and  presented  a  fig¬ 
ure  showing  the  oceanic  and  continental  shear 
velocity  distributions  derived  from  calculations 
for  flat  layers.  These  authors  emphasized  that  a 
proper  treatment  of  the  problem  must  include 
sphericity. 

Theoretical  studies  have  recently  been  made 
which  include  allowance  for  sphericity,  gravity, 
gradients,  and  a  liquid  core.  The  spheroidal  os¬ 
cillations  have  been  treated  theoretically  by  Al- 
terman  ct  al.  [1961],  Bolt  and  Dorman  [1961], 
and  others.  The  torsional  oscillations  have  been 
investigated  by  Jobert  [1959,  1960a,  6],  Sato 
et  al.  [1960a,  6],  Pekeris  et  al.  [1961],  MacDoti- 
ald  and  Ness  [1961],  and  others.  Jobert  [19606] 
presented  Love  wave  phase  and  group  velocity 
cur\'es  for  a  modified  Gutenberg  mantle  which 
was  considered  appropriate  for  oceanic  areas. 
This  calculation  for  a  spherical  earth  produced 
a  fairly  flat  group  velocity  curve  with  a  majd- 
muin  of  approximately  4,5  kni/sec  at  30  seconds 
and  a  minimum  of  aj>proximatcly  4.4  km/sec 
near  3  minutes.  Although  no  comp.ariEon  was 
m.ade  with  observed  data,  the  maximum  .and 
the  minimum  lie  slightly  outside  the  range  of 
ob-served  group  velocities  for  the  fundamental 
oceanic  Love  mode. 

Some  of  the  theoretical  studies  cited  abovc- 
were  .stimulated  by  the  occurrence  of  the  great 
Chilean  ear'thquoke  of  May  22,  1930,  which 
excited  the  free  vibrations  reported  by  A /sop 
et  al.  [1961],  BeniofJ  et  al.  [1961],  Bogert 
[1961],  Ness  et  al.  [1961],  Brune,  Benioff,  and 
Ewing  [1961],  and  others.  These  studies  of  the 
Chilean  earthquake  prod'-'-ed  new  and  precise 
dctermin.ations  of  the  dispersion  of  long-period 
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Love  and  uaive.*'.  For  periods  greater 

(liaji  several  Imndrod  Hecoiidh',  the  Ktecp  rise  in 
])hasc  and  group  ve!oeitir,«,  whirh  had  been  pre¬ 
dicted  by  theory,  was  eonrirmed  by  observations. 
This  steel)  ri.‘-e  is  paused  by  the  penetration  of 
the  wave  disturbance  info  tlie  lower  half  of  the 
mantle. 

Com'fndntionnl  proccdHrvK.  I<ove  wave  dis- 
j)ersion  for  a  spherical  {'arth,  consisting  of  a 
large  number  of  splu'fical  shells,  has  been  com¬ 
puted  using  an  IHM  7()!H)  version  of  a  nu¬ 
merical  iteration  program  desj'riiied  by  Sain 
ct  al.  [idOOo,  6].  This  firograin  lia.s  been  u.sed 
to  calculate  free  iwriods  of  torsional  oscil¬ 
lation,  uhaso  velocity,  group  velocity,  particle 
amplitude,  and  .strain  di.slribution  as  a  func¬ 
tion  of  depth  for  each  of  the  azimuthal  mode 
numbers  n.  Group  velocity,  which  is  calculated 
by  the  polynomial  method,  is  about  two  orders 
of  magnitude  more  preei.se  than  the  observa- 
tion.al  data.  With  this  program  dispersion  can 
bo  computed  for  a  spherical  or  a  flat-layered 
model;  the  same  physical  parameters  and  nu¬ 
merical  procedures  can  be  used  in  both  cases. 
These  oalculat!''ns  for  flat  layers  were  compared 
with  those  made  with  the  aid  of  the  Thomson- 
Haskcll  [Haskell,  lOliS]  matrix  program,  PV7, 
w’ritten  by  Dorman  [Oliver  and  Dorman,  1981], 
and  the  results  agree  at  all  periods  to  better 
than  0,0001  km/scc.  In  addition,  the  matrix 
method  was  combined  with  the  earth-flattening 
approximation  [Alterman  rt  al.,  1981]  as  a  fur¬ 
ther  test  of  the  caleulntion.i  for  the  spherical 
model.  All  methods  give  very  similar  funda¬ 
mental  Love  mode  phase  velocities  for  periods 
shorter  than  8  seconds.  For  period.s  shorter  than 
50  seconds  the  earth-flattening  approximation 
ngrees  to  bettor  than  0.01  kni/.«ec  with  the  fun¬ 
damental  Love  mode  plin.se  .and  group  velocities 
determined  by  tlie  ealeuhifion.s  for  a  spheric.'d 
earth. 

A  jirogram  written  by  L.  E.  Alsop  was  n.scd 
in  romputing  oceanic  llaj'Ieigli  wave  free  pe¬ 
riods,  7)lia.se  and  group  veloeitie.':,  particle  am- 
?(';t,'.ide«,  -and  gr.avitational  di.sturhnnce  for  e.aeh 
a;-.imut!ial  mode  of  a  gravitating,  .spherical  earth 
with  a  liquid  surface  layer.  This  routine  is  simi¬ 
lar  to  one  dc.scribed  by  Alsnp  [1902].  The 
oreatiic  Tv.a.yleigh  mode  phase  and  group  veloci¬ 
ties  which  ro.sulted  were  tlien  u.sed  to  check  the 
correction  for  nnhericity  for  ranlinrntal  Tfay- 
It'igli  w.'ives'  reported  by  Balt  end  Dorman 
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[1981],  which  wa.s  applied  to  Dorman’s  oceanic 
ease  8999  by  Kno  ct  al.  [1962]. 

For  periotk  between  .90  and  150  seconds,  tlie 
phase  and  group  velocities  calculated  for  c.a.se 
122  agree  with  the  theoretical  values  for  ca.se 
S099S  (spherical)  of  Kuo  et  al.  to  within  0.01 
km/scc.  Rayleigh  wave  ea]cuIation,s  for  a  flat- 
l.ayered  version  of  case  122  also  agree  with  Dor¬ 
man’s  published  values  for  case  8099  to  within 
0.01  km/sec  over  the  same  range  of  periods 
For  periods  shorter  than  50  seconds  the  eartli- 
rlatteiiirig  approximation  for  case  122  agrees  to 
better  than  0.01  km/sec  with  the  Rayleigh  mode 
phase  and  group  velocities  for  a  sphorical  earth. 

A  matrix  ptograni  for  liquid  and  solid 
layers  was  used  in  computing  Rayleigh  w'avc 
particle  amplitudes.  The  method  is  similar  to 
that  reported  by  Dorman  and  Prentiss  [1960] 
and  Dorman  [1962]. 

Oceanic  case  122  cermpared  with  other  studies 
of  the  crust  and  mantle.  For  the  computations 
of  spherical  Love  and  Rayleigh  waves  for  the 
mode!  reported  in  this  study,  the  velocities  and 
densities  for  the  crust  and  mantle  were  chosen 
to  be  in  accord  with  recent  geophysical  studies 
of  oceanic  area.s.  The  velocities  for  the  sedimen¬ 
tary  and  crustal  layers  of  case  122  were  indi¬ 
cated  by  the  recent  work  of  Oliver  and  Dorman 
[1961]  and  are  consistent  with  seismic  refrac¬ 
tion  work  at  sea  as  reported  by  Raltt  [1956] 
and  Ewing  and  Ewing  [1959].  Crustal  densities 
adopted  for  case  122  conform  to  the  Nafo  and 
Drake  em.piricn.1  velocity-density  curve  reported 
by  Talwani  et  al.  [1959].  The  layer  parameters 
for  the  homogeneous  isotropic  layers  which  con¬ 
stitute  ca,se  122  are  given  in  Table  1.  For  the 
crust,  these  layers  correspond  closely  in  both 
thickness  and  velocity  to  those  indicated  by 
seismic  refraction.  Ewing  and  Ewing  [1961], 
using  a  seismograph  on  the  ocean  floor,  obtained 
sub.ccrn.'dal  shear  velocity  mea.sv.rements  at  sm.vil 
distances.  These  autliors  consider  tlieir  me.tsure- 
ments  to  be  appropriate  only  for  the  upper  sur¬ 
face  of  the  rn.antle,  and  for  thi.s  rca.soii  their 
value  of  mantle  .shear  velocity  w.os  not  u  cd  for 
case  122. 

For  Lie  mantle  the  layer  parameter,  for  ca-;- 
122,  pre.aented  in  Table  1,  are  identical  to  tliov' 
for  oceanic  c.ase  8099  cf  Dorman  et  .al.  The 
homogeneous  l.a.yors  for  the  mantle  were  rc- 
uuired  bv  Dorman  et  al.  Iiecau,-.;'  of  Iheir  o’cthod 
of  caleuiatioi).  'Dic.se  layer:-'  .••rc  a  .-depped  ap.- 
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TABLE  1 .  Physical  Parameters  for  the  Homogeneous,  Isotropic,  Spherical  Shells  of  Case  122 
Thickness  (^),  compressional  velocity  a  (km/sec),  shear  velocity  /?  (km/sec),  and  density 
-  ®370  km  and  core  radius  =  3470  km  for  the  calculations  for  the 
spher  cal  model.  Deepest  layer  represents  the  half-space  for  calculations  for  the  flat-layered 


LAYER 

NO, 


THICKNESS 

KM 

5.0 

0.5 

1.5 

,5.0 

49.0 

i6o.o 

ico.o 

90.0 

90.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

lOOoO 

100.0 

100.0 

100.0 

100.0 

100.0 


KM/S 

1.52 

2.00 

4.80 

6.90 

8.49 

8.81 
9.112 

.lil 

10.85 
11.12 
11.33 

11.49 
11.64 
11,78 

11.92 

12.06 

12.19 

12.^3 

12.46 

12,59 

12.72 

12.85 

12,97 

13.09 

13.21 

J3.33 

13,40 

13.^3 

13.61 

ii64 


p 

KM/S 

0.0000 

0.5000 

2.7700 

3.9800 

4.6125 

4.3000 
4.6000 
4.8000 
5.1925 
5.4925 
5.7900 
6.0300 
6.2000 
6.3150 

6.4000 

6.46SO 

6.5310 

6.5910 

6.6490 

6.7o4o 

6.7550 

6.8050 

6.8525 

6.8975 

6.9%o 
6.9950 
7.0450 
7.0950 
7. 1425 

7.1875 

7.2325 

7.2775 

7.3000 


p 

GM/CC 

1.0300 

1.9000 
2.5400 

2.9000 
3.3,400 
3.4425 
3-5265 

3,6o4o 

3.7650 

4.0100 

4.2300 

4.4100 

4.5450 
4.6400 
4.7100 
4.7700 
4.8275 
4.8825 
4.9400 
5.0000 
5.0550 
5. 1050 

5.1575 

5.2075 

5.2650 

5.3150 


5.4650 

5.5150 

5.5625 

5.6075 

5.6550 


proximation  tp  volocu  ics  and  donsU!’.- :  v/iucji  .-ire 
l)a.scd  on  tho.S'C  ori^inaHy  detcnnini’d  *>001  (-'riii- 
quake  body  wave  studies.  In  the  present  study, 
the  homogeneous  i.sotropic  layers  of  case  8099 
of  Dorman  et  a.!,  have  been  retained  for  case 
122. 

Ob.servation.s  of  the  travel  times  of  P  and  .8 
from  nuclear  explo.sions  in  the  Pacific  by  Carder 
and  Bailey  [1958],  Kogan  [1960],  and  others 
have  confirmed  the  telc.seism.ic  travel  time  curves 
remarkably  well.  For  source.s  in  the  Pacific  and 
recording  sial  ions  on  continents,  P  arrives  2  sec¬ 
onds  early  and  8  is  recorded  4  -|-  seconds  late, 
as  compared  witli  the  Jeff reys-Bullen  tables.  Ex¬ 


cept  ior  one  observation  of  S  near  17'’  [Pomeroy, 
1902  j,  no  i  ravol-time  data  are  available  from  nu¬ 
clear  explosions  in  oceanic  areas  at  distances 
less  thiin  34"  for  S  and  3]4°  for  P.  The  short- 
{Mjried  S  wave  observation  near  17°  indicates  a 
minimum  upper-mantle  shear  velocity  of  4.7 
km/sec  beneath  the  ocean  basin  between  Eni- 
wetok  and  Guam.  No  oceanic  refraction  re¬ 
sults  exist  for  distances  greater  than  1°.  These 
precise  results  must  be  further  amplified  in 
order  that  they  may  gi\e  information  about 
the  suboceanic  upper  mantle. 

It  is  instructive  to  compare  the  mantle  ve¬ 
locities  determined  from  surface  wave  studies 


MANTLE  SHEAR  WAVE  VELOCITIES 


of  oceanic  and  continental  areas.  The  sub- 
oceanic  mantle  distributions  to  be  considered 
are  case  122,  case  8099  of  Dorman  et  al.  [1960], 
and  case  6EGHP1  of  Aki  and  Press  [196ll 
These  distributions  are  very  similar;  all  contain 
a  channel  with  a  minimum  shear  velocity  of 
4.30  km/sec.  As  noted  above,  cases  122  and 
j  8099  are  identical  in  the  mantle.  The  upper  and 
lower  channel  boundaries  of  case  6EGHP1  are 
not  as  sharp  as  those  for  cases  122  and  8099, 
since  case  8EGHP1  was  obtained  by  a  modifica- 
tion  of  a  Gutenberg  instead  of  a  Lehmann  con¬ 
tinental  distribution. 

For  the  continents  the  Gutenberg-Bullen  A 
model,  case  6EGH  of  Aki  and  Press  [1961],  and 
the  distribution  CANSD  determined  from  the 
dispersion  of  Love  and  Rayleigh  waves  across 
the  Canadian  shield  [Brune  and  Dorman,  1962] 
will  now  be  conadered.  The  Gutenberg-Bullen 
A  model  has  been  used  to  explain  the  observed 
mantle  Love  and  Ra3deigh  wave  data  from  the 
great  Chilean  earthquake  and  the  southeast 
Alaska  shock  of  July  10,  1958  (e.g.,  Brune, 
BeniofJ,  and  Ewing  [1961]  :  Pekeris  et  al.  [1961], 
r  c.).  Case  6EGH  is  a  model  based  on  one  of 
Ciitcnberg’s  idter  velocity  distributions.  The 
principal  difference  between  case  6EGII  and  the 
Gutenberg-Bullen  A  model  occurs  between  50 
and  100  km,  where  case  BEGE  has  lower  shear 
velocities.  The  Canadian  shield  case,  CANSD, 
is  characterized  by  higher  velocities  than  either 
of  these  models.  The  shear  velocity  from  the 
Mohorovicic  discontinuity  to  a  depth  of  115  km 
is  4.72  km/sec.  The  low- velocity  channel,  which 
is  200-km  thick,  consists  of  a  4.54-km/Bec  layer 
overlying  a  4.51-km/sec  layer  of  equal  thick¬ 
ness. 

In  general,  surface  wave  studies  to  date  have 
shown  that  the  channel  velocity  is  lower  be- 
•  neath  the  oceans  than  beneath  the  continents. 
In  addition,  with  the  exception  of  the  study 
of  Aki  and  Press  [1961],  several  of  these  inves¬ 
tigations  have  indicated  that  the  oceanic  chiiii- 
nel  velocities  occur  nearer  the  surface.  Below 
the  M  di.scontinuity,  the  physical  parameters  for 
the  continental  model  6EGH  and  the  oceanic 
model  6EGHP1  of  Aki  and  Press  differ  only  in 
the  value  of  shear  velocity  chosen  for  the  low- 
velocity  zone.  Model  6EGH  was  compared  with 
only  one  seismogram  for  periods  near  tlie  con¬ 
tinental  Air>'-phase  maximum.  Further  study 
and  a  comparison  of  various  models  with  more 
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continental  data  arc  needed  in  order  to  de¬ 
termine  whether  channel  velocities  usually  oc¬ 
cur  at  greater  depths  beneath  the  continents 
than  beneath  the  oceans. 

Love  wave  phase  and  group  velocities.  Ex¬ 
perimental  and  theoretical  fundamental  mode 
phase  velocities  for  oceanic  Love  waves  of  pe¬ 
riods  up  to  800  seconds  are  presented  in  Rgure  1. 
Experimental  data  appropriate  for  oceanic  areas 
have  been  taken  from  Satd  [1958],  as  corrected 
for  the  polar  phase  shift  [Brune,  Nafe,  and  Al- 
sop,  1981;  Brune,  Beniojf,  and  Ewing,  1961; 
Benioff  et  al.,  1961 ;  MacDonald  and  Ness,  1961] . 
Of  these,  Sato's  data  have  the  greatest  per¬ 
centage  of  oceanic  path  and  hence  his  data  are 
considered  more  reliable.  Four  theoretical  curves 
are  shown;  one  is  the  result  of  a  calculation  for 
flat  layers,  and  the  other  three  are  for  a  spherical 
earth.  Numerical  results  of  the  calculations  of 
Love  waves  for  the  spherical  version  of  case 
122  appear  in  Table  2.  For  azimuthal  mode  num¬ 
bers  less  than  360,  the  second  and  third  layers 
(shear  velocity  =  0.5  and  2.77  km/sec,  respec¬ 
tively)  were  deleted  and  replaced  by  a  single  2- 
km  layer  with  a  shear  velocity  of  2.77  km/sec 
and  a  density  of  2.54  g/cm*.  A  detailed  com¬ 
parison  of  these  computations  with  data  will  be 
made  in  conjunction  with  the  material  presented 
in  Figure  2. 

Calculated  group  velocities  for  two  flat-layered 
and  two  spherical  models  are  compared  with 
group  velocity  data  appropriate  for  oceanic 
areas  in  Rgure  2.  The  sources  of  these  data  are 
indicated  in  the  figure.  In  the  long-period  Chile 
to  Nana  data,  those  represented  by  stars  are  the 
more  reliable,  according  to  Brune,  Beniojf,  and 
Ewing  [1961] ,  because  the  effects  of  initial  phase 
at  the  source  have  been  removed.  The  symbol  G 
represents  values  which  were  picket!  from  the 
summary  curve  of  I^ove  wave  data  for  the  Pa¬ 
cific  in  Figijre  1  in  the  third  paper  of  Gutenberg 
and  Richter  [1936,  p.  89].  The  data  indicated 
by  squares  are  from  a  southeast  Pacific  shock 
which  occurred  on  November  14,  1958,  at  05h 
04m  253  GCT  at  36.0''S,  102.8*W,  a.s  recorded 
at  En’lett  otaiicn,  Antarctica. 

Several  observations  may  be  made  from  a 
studj'  of  Figures  1  and  2.  The  calcvilation  of  case 
122  for  flat  layers  departs  significantly  from 
the  Love  wave  phase  end  j'roup  velocity  data 
in  Figures  1  and  2.  In  contrast,  when  the  com¬ 
putation  is  done  foi  a  sphericn!  earth,  the  phase 
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Fig.  1.  Phase  velocities  for  fundamental-mode  oceanic  Love  waves.  Calculation  for  a 
spherical  earth  (solid  line)  agrees  more  closely  with  observational  data  than  calculaticm  for 
a  flat-layered  model  (dashed  line) .  Parameters  for  case  122  in  Table  1. 


velocities  .''^ree  with  a'^iiiable  data  for  periods 
between  60  and  800  seconds  to  within  a  few 
hundredths  of  a  kilometer  per  second.  The  cor¬ 
responding  group  velocities  compare  well  with 
existing  data  over  an  even  broader  range  of 
peiiods.  The  agreement  of  theoretical  wave  ve- 
iocitir-s  for  case  122  with  observational  data  over 
such  a  broad  range  of  [leriods  indicates  that  the 
physical  parajueters  chosen  are  close  to  the  ac¬ 
tual  ones  in  the  earth  to  depths  of  at  least  sev¬ 
eral  hundred  kilometers  beneath  the  ocean. 
These  results  may  be  further  improved  by 
minor  modifications  of  the  model.  New  and  more 
precise  data  are  also  needed. 

If  no  calculations  for  a  spherical  earth  were 
available,  one  might  tend  to  modify  the  physical 
parameters  of  case  122  in  order  to  explain  the 
Love  wave  phase  and  group  velocity  data.  One 
of  the  better  Love  wave  caicnlations  for  flat 
layers  is  the  c-ase  marked  Oceans  VTII  [Larzrft’s- 
mtLP.  ct  ni.  1059].  The  shear  velocity  distribu¬ 
tion  for  this  case  was  presented  by  Eiving  and 
Landitivian  [1061].  The  calcvilation  of  phase  and 
group  velocities  for  this  model  was  the  first 
f|uanlitative  attempt  to  explain  the  disjxTvsion 
of  oceani;-  Love  waves  with  a  mathematicai 
model  eontnining  a  low-velocity  ehannei.  This 


model  indicated  that  oceanic  charnel  velocities 
are  smaller  and  come  nearer  to  the  surface  than 
continental  channel  velocities.  As  was  discussed 
above,  several  recent  studies  corroborate  this 
result.  In  comparison  with  case.s  122  and  8099, 
the  shear  velocity  for  Oceans  VIII  is  higher  by 
several  tenths  of  a  kilometer  per  .second  in  the 
lower  part  of  the  channel.  The  differences  be¬ 
tween  the  shear  velocity  distributions  for  Oceans 
VIII,  derived  from  Love  wave  dispersion,  and 
case  8099  of  Dorman  et  al.,  determined  from 
Rayleigh  wave  data,  were  one  of  the  discrep¬ 
ancies  noted  by  Ewing  and  Landisman.  Similar 
discrepancies  between  Love  and  Ralyeigh  wave 
results  have  been  the  basis  of  recent  speculations 
concerning  anisotropy  in  the  upper  mantle  [An¬ 
derson,  1961,  1962;  Anderson  and  Harkrider, 
1962].  Since  case  122  agrees  with  the  data 
when  calculation  for  a  spherical  earth  is  per¬ 
formed,  no  recourse  to  a  difference  in  SH  .-.nd 
5V  velocities,  i.e..  snlstvopy,  reouired.  Thus 
one  di.stribution  of  homogenecus  isotropic  lay¬ 
ers  can  satisfy  both  T.ove  and  Rayleigh  wave 
data. 

In  I'bgurrs  1  and  2,  the  higher  velocitier,  at 
longer  j>t‘ricds  have  been  iricjisured  using  data 
from  the  great  Chilean  earrhquake.  T'.e  t  iicoreti- 
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Fig.  2.  Group  velocities  for  fundamental-mode  oceanic  Love  waves.  Calculation  for 
sphf'~icn,l  earth  case  122,  compared  with  observational  d.afa  from  10  to  SOO  seconds  period.  For 
diaou.ssion  see  te”*. 


cal  curves  for  the  Jeffreys-Bullcn  and  Guten- 
berg-Bullen  A  continental  models  shown  in  Fig¬ 
ures  1  and  2  agree  quite  wcl!  with  the  results 
of  case  122  at  long  jjeriods. 

For  periods  shorter  than  12  seconds  the  lower- 
velocity  oceanic  crustal  materials,  in  particular 
the  sediments,  produce  a  sharp  drop  in  group 
velocity  as  the  period  decreases.  For  periods  be¬ 
tween  10  and  approxim.ately  50  seconds,  small 
segments  of  continental  path  can  lower  the  ob- 
■served  grouj)  velocities.  For  the  data  fso::'!  Ilai- 
lett  station,  a  correction  was  appliet:  for  the  G 
per  cent  of  the  path  that  lies  within  the  continent 
of  Antarctica. 

The  wave  guide  of  the  G  wave.  Tlie  well- 
known  pulselike  character  of  the  G  wave,  v.-hich 
i.-;  best  observed  for  oceanic  paths,  is  produced 
by  the  nearly  constant  group  velocity  of  the 
fundamental  oceanic  Love  mode  for  periods  as 
great  as  300  seconds,  as  shown  in  Figure  2. 
In  Ills  study  of  the  wave  guide  of  the  G  wave, 
Sntd  [1958]  said  that  ‘this  phase  corresponds  to 
a  Love  wave  traveling  in  the  mantle  and  con- 
trollffl  by  th<‘  gnid'cnt  in  ‘-I'car  velocity’  between 


depths  of  100  and  several  hundred  kilometers. 

Wc  may  further  study  the  wave  guide  of  the 
G  wave  with  the  aid  of  the  Love  wave  particle 
amplitudes  computed  for  a  .spherical  earth  (Fig¬ 
ure  3).  For  reference,  the  distribution  of  shear 
velocity  with  depth  is  shown  at  the  left.  The 
particle  amplitudes  are  normalized  to  unity  at 
the  free  surface.  The  particle  amplitude  for  T 
~  59.07  seconds  is  representative  of  particle 
amplitudes  for  the  G  wave  for  periods  between 
GO  and  300  seconds.  For  the.se  periods  the  free 
.surface  and  the  .steep  gradient  in  shear  velocity 
down  to  approximately  900  km  form  the  wave 
guide  of  the  G  wave.  For  periods  greater  th-an 
300  seconds,  the  particle  motion  penetrates  the 
high-speed  lower  mantle,  and  observed  Love 
wav^o  pnaGS  and  group  vclocitie.''  incre.nse  sharply, 
.03  had  been  predicted  by  thcorjs  Similarly,  for 
periods  grealei  lh..ri  250  seconds,  ob, served  and 
theoretical  Rayleigh  wave  ]>liase  and  group  ve¬ 
locities  rise  even  more  abniptly  than  those  for 
Love  waves.  For  these  wavelengths,  Rayleigh 
waves  penetrate  more  deeply  into  the  earth. 

In  Figures  1  and  2,  the  i>haso  and  group  vc- 
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Fig.  3.  Particle  amplitude  profiles  as  a  function  of  depth  for  the  fundamental  Lo'/e  mode 
of  case  122,  calculation  for  spherical  model.  Amplitudes  normalized  to  unity  at  the  free  surface. 
Secondary  maximum  develops  when  phase  velocity  C  begins  to  exceed  4.3  km/sec,  the  channel 
velocity. 


locities  for  case  122,  as  calculated  for  a  spherical 
earth  and  for  fiat  layer.s,  remain  approximately 
parallel  for  periods  between  12  and  150  seconds. 
For  periods  shorter  than  10  seconds,  the  two 
phase  velocity  curves  converge;  for  periods 
*  shorter  than  8  seconds,  they  are  indistinguish¬ 
able.  As  discussed  in  the  section  on  computa¬ 
tional  procedures,  both  numerical  iteration  and 
^  the  matrix  method  were  used  for  a"  fhese  cal¬ 
culations  and  were  found  to  be  in  good  agree¬ 
ment. 

Within  the  period  range  where  the  phase  and 
group  velocity  curves  are  parallel,  the  difference 
between  the  calculations  for  the  spherical  earth 
and  for  the  flat-layered  versions  of  case  122  is 
approximately  0.8  km/sec.  In  simplest  terms, 
this  increase  in  velocity  for  the  spherical  earth 
can  bo  attributed  to  a  shorter  path.  For  the 
spherical  earth,  if  we  imagine  that  the  tvaves 
travel  predominantly  in  the  low-velocity  chan¬ 
nel  iastcad  of  along  the  surface,  the  path  is 


shortened  by  a  factor  approximately  equal  to 
the  ratio  of  the  effective  channel  depth  to  the 
ear  ill's  radius  (M.  Ewing,  personal  communica- 
cation).  Thus  the  velocity  of  the  waves  is  in¬ 
creased  by  approximately  (s/a)  •  where  s 
is  the  effective  depth  to  the  channel,  a  the 
radius  of  the  earth,  and  V,  the  velocity  in  the 
channel.  In  the  present  case,  with  a  channel  ve¬ 
locity  of  4.30  km/sec,  the  difference  of  0.08  km/ 
sec  in  Figures  1  and  2  indicates  an  effoctive  chan¬ 
nel  depth  of  120  km.  This  depth  lies  well  within 
the  channel  for  case  122.  With  this  difference  of 
path  length  in  mind,  it  is  easy  to  understand 
hov/  a  spherical  model  with  a  channel  velocity 
of  4.3  km/sec  devoioiis  the  observed  G  wave 
velocity  of  approximately  4.4  km/sec.  Press  and 
Exving  [1956]  and  Press  [1959]  have  suggested, 
without  detailed  calculation,  that  the  compari¬ 
son  of  the  4.7-km/sec  mantle  shear  velocity,  de¬ 
rived  from  .seismic  refraction  nieasuremeats, 
with  the  4.4-kin/scc  G-wave  velocity  indicates 
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the  presence  of  a  low-velocity  channel  in  the 
upper  mantle. 

To  understand  the  reasons  for  the  abrupt  di¬ 
vergence  of  the  phase  velocity  curves  for  spheri¬ 
cal  and  flat  models  near  10  seconds,  we  may 
consider  the  diagrams  of  Love  wave  particle  am¬ 
plitude  for  a  spherical  earth  shown  in  Figure  3. 
At  the  shortest  period,  8.03  seconds,  the  motion 
is  confined  to  the  crust  and  to  the  high-speed 
layer  just  below  the  crust.  In  this  part  of  the 
spectrum,  the  calculations  for  spherical  and  flat 
models  give  nearly  identical  phase  and  group 
velocities  and  particle  amplitudes.  When  the  pe¬ 
riod  becomes  greater  than  8  seconds,  the  phase 
velocity  increases  and  the  waves  penetrate 
farther  into  the  high-speed  portion  of  the  man¬ 
tle.  At  the  period  where  the  phase  velocity  be¬ 
comes  equal  to  the  channel  velocity,  drastic 
changes  in  particle  amplitude  begin  to  occur. 
Within  a  period  range  of  approximately  1  sec¬ 
ond,  the  wave  which  was  confined  above  the 
channel  emerges  as  a  disturbance  trapped  in 
the  entire  low-velocity  zone.  This  penetration 
of  the  surface  waves  through  the  high-speed 
mantle  into  the  low-velocity  channel  is  analogous 
to  many  problems  of  barrier  penetration  in 
quantum  mechanics. 

Near  10  seconds,  where  the  fundamental 
oceanic  Love  mode  penetrates  the  chann&l,  the 
particle  amplitudes  for  the  spherical  a;nd  flat¬ 
layered  versions  of  model  122  differ  and  the 
phase  velocity  curves  diverge.  Calculations  of 
particle  amplitudes  for  the  flat-layered  version 
show  more  severe  trapping.  The  lower  phase  ve¬ 
locity  and  larger  channel  amplitudes  for  this 
case  may  be  attributed  to  a  more  effective  low- 
velocity  channel. 

The  particle  amplitude  curves  for  the  spheri¬ 
cal  earth  version  of  case  122  at  T  =  10.34  and 
15.23  seconds  in  Figure  3  illustrate  trapping  and 
also  exhibit  secondary  maximums.  For  case  122, 
the  largest  secondary  maximum  occurs  near  15 
seconds  when  the  calculation  is  performed  for  a 
spherical  earth.  At  these  periods,  the  phase  and 
group  velocities  and  the  particle  amplitudes 
which  result  from  the  earth-flattening  approxi¬ 
mation  are  extremely  similar  to  those  calculated 
for  the  spherical  earth.  Thus,  for  the  spherical 
earth,  the  simplified,  explanation  discussed  pre¬ 
viously  for  the  relation  between  the  channel 
velocity  and  the  G  wave  velocity  is  justified  by 
the  study  of  particle  amplitude  trapping  and 


secondary  marimums  in  the  low-velocity  chan¬ 
nel  for  case  122; 

There  are  two  requirements  for  the  develop¬ 
ment  of  secondary  maximums  in  Love  wave 
particle  amplitudes.  First,  the  phase  velocity 
must  be  slightly  higher  than  the  channel  ve¬ 
locity  but  less  than  the  velocity  of  the  higher- 
speed  regions  adjacent  to  the  channel.  Second, 
the  vertical  wavelength  must  be  less  than  or 
comparable  to  the  channel  thickness. 

In  a  study  of  oceanic  Love  waves,  Johert 
[19606]  showed  that,  for  n  =  600,  the  fundar 
mental  mode  particle  amplitude  develops  a  sec¬ 
ondary  maximum.  Said  et  al.  [19606]  also  ob¬ 
served  that,  for  models  with  a  low-velocity  layer, 
the  higher-mode  particle  amplitudes  can  be  much 
greater  in  the  channel  than  at  the  surface. 

It  is  evident  that  secondary  maximums  will 
play  an  important  role  in  the  problem  of  nor¬ 
mal  mode  excitation  by  earthquakes  at  depth. 
Ccdoi  [1953,  1954]  and  Gutenberg  [1954, 1955] 
reported  that  earthquakes  within  the  low-veloc¬ 
ity  channel  tend  to  excite  the  phases  P«  and 
S«  quite  easily.  Bolt  and  Dorman  [1961]  calcu¬ 
lated  the  group  velocity  of  the  first  shear  mode 
for  a  Gutenberg  velocity  distribution.  They 
found  a  maximum  in  the  group  velocity  curve 
at  4.54  km/sec  for  a  period  of  25  seconds  which 
agrees  well  with  the  period  and  velocity  of  the 
phase  S.  described  by  Caloi  and  Gutenberg  and 
the  phase  Sn  of  Press  and  JEwing  [1955]. 

The  fundamental  oceanic  Love  mode  for  case 
122  has  a  secondary  noaximum  of  particle 
amplitude  at  short  periods.  For  the  Canadian 
shield,  Brune  and  Dorman  [1962]  have  presented 
a  case,  CANSD,  for  which  the  fundamental  Love 
mode  is  not  associated  with  channel  waves.  A 
secondary  maximum  does  not  develop  in  this 
case,  nor  in  the  Gutenberg-Bullen  A.  model,  since 
the  phase  velocity  does  not  exceed  the  lowest 
channel  velocity  until  the  vertical  wavelength 
is  greater  than  the  channel  thickness.  For  these 
models,  however,  the  first  shear  and  the  secpnd 
Love  modes  are  much  like  the  fundamental  Love 
mode  for  the  pceans  in  that  a  secondary  maxi¬ 
mum,  which  always  characterizes  these  modes, 
is  confined  to  the  low-velocity  zone.  All  these 
modes  are  associated  with  channel  waves.  The 
higher  modes  offer  a  means  of  refining  the  knowl¬ 
edge  of  the  mantle  beneatli  various  regions. 

The  oceanic  Rayleigh  mode.  Oceanic  Ray¬ 
leigh  mode  phase  and  group  velocities  as  ob- 
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served  for  Pacific  paths  for  periods  between  20  a  program  written  by  L.  E.  Alsdp  of  the  Lamdnt 

and  150  seconds  have  been  presented  by  Kuo  Geological  Observatory  was  used  to  check  this 

et  oZ.  [1962],  AH  the  data  presented  in  Figures  27  correction.  The  exact  calculation  for  case  122 

and  28  of  their  study,  with  the  exception  of  those  .agrees  with  curves  8099S  to  within  OiOl  km/sec. 

labeled  P13  and  P14,  may  be  considered  repre-  Thus,  the  observed  oceanic  Rayleigh  mode  phase 

sentative  of  the  dispersion  for  Pacific  Ocean  and  group  yelpcities  for  Pacific  paths  are  in  close 

basins.  The  paths  indicated  by  data  points  P13  agreement  with  those  calculated  for  case  122.  In 

and  P14  cross  the  Melanesia-New  Zealand  re-  the  preceding  section  it  was  shown  that  this 

gion  of  the  southwest  Pacific  and  hence  will  not  case  agrees  with  dispersion  data  for  oceanic 

be  included  in  the  present  discussion.  The  scat-  Love  waves.  These  i-nn  pnrisons  indicate  that 

ter  of  the  data  below  30  seconds  can  be  at-  case  122  must  not  dr  iKirt  greatly  from  the  condi-. 
tributed  to  differences  in  water  depth.  tions  to  be  foui  1  in  tlio  niantle  beneath  the 

The  observed  phase  and  group  velocities  of  oceans,  since  it  satisfies  both  Love  and  Ray- 
Kuo  et  al.  were  compared  with  theoretical  leigh  wave  data. 

curves  designated  8099S.  These  curves  were  Rayleigh  wave  amplitudes  for  case  122,  corn- 
obtained  by  correcting  the  phase  velocities  for  parable  with  the  amplitudes  for  Love  waves  in 

case  8099  of  Dorman  et  al.  [J960]  according  to  Figure  3,  are  shown  in  Figure  4.  The  distribution 

the  spherical  correction  formula  for  the  Ray-  of  shear  vele-city  with  depth  is  again  pveii;  at  the 

leigh  wave  phase  velocities  reported  hy  Bolt  and  left  side  of  the  figure.  The  vertical  displace- 

Dorman  [1961]  for  the  Gutenberg-Bullen  A  ments  are  normalized  to  unity  at  the  water  sur- 

continental  model.  face.  Particle  amplitudes  for  case  122  are  shown 

As  was  noted  in  the  section  on  computations,  for  one  calculation  for  a  spherical  earth  and  three 

SHEAR  VELOCITY  PARTICLE  AMPLITUDES 

(KM /SEC) 


Fig.  4.  Particle  amplitude  profiles  as  a  function  of  depth  for  the  Rayleigh  mode  of  case 
122.  Vertical  displacement  V,  normalized  to  unity  at  the  free  surface.  H  is  horizontal  displace¬ 
ment,  H  and  V  of  same  sign  indicate  retrograde  motion.  For  a  given  wayelength,  X,  particle 
amplitudes  are  the  same  for  calculations  for  flat  and  spherical  models.  To  first  order,  the  group 
velocities  are  also  the  same. 
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calculations  for  flat  layers.  No  secondary  maxi¬ 
mum  or  channel  wave  develops  in  the  low-veloc¬ 
ity  channel  for  the  oceanic  Rayleigh  mode  of 
case  122. 

When  calculations  of  Rayleigh  mode  particle 
amplitudes  for  spherical  and  flat-layered  versions 
of  case  122  are  compared  for  a  given  wavelength 
(Figure  4),  it  is  found  that  the  particle  ampli¬ 
tudes  are  nearly  the  same.  For  periods  shorter 
than  100  seconds,  the  group  velocities  agree  to 
within  0.02  km/sec.  Bolt  and  Dorman  [1961] 
also  reported  that  the  Rayleigh  mode  group  ve¬ 
locities  calculated  for  spherical  arid  flat-layered 
versions  of  their  models  are  nearly  the  same, 
even  though  the  phase  velocities  are  different. 
To  extend  the  discussion  of  the  relation  between 
phase  and  group  velocities  as  a  function  of 
wavelength  for  a  spherical  and  a  flat-layered 
medium,  let  us  consider  one  of  the  well-known 
expressions  which  relates  the  group  velocity  U 
to  the  phase  velocity  C  for  a  certain  wave¬ 
length  A: 

U.-=  G.-\  dC./d\  (la) 

where  the  subscript  3  indicates  a  spherical  earth. 
Similarly,  for  a  flat  earth, 

Uf=C,-\  dCf/dX  (16) 

If  the  group  velocities  and  wavelengths  are  ap¬ 
proximately  equal,  we  may  subtract  (16)  from 
(la)  and  obtain 

0  (C.  -  €f)  -  X[dC./dA  -  dC,/dK]  (2) 

This  may  be  rearranged  as 

(C.  -  Cf)/\  ^  d{C.  -  C,)/d\  TH  m  (3) 

where  m  is  a  constant.  Alternatively,  we  may 
write 

C.  «  C,[l  -I-  T  d{G.  -  G,)/dX] 

Gf[l  -f  mT] 

The  last  of  equations  4  is  identical  to  the  em¬ 
pirical  formula  (12)  of  Bolt  and  Dorman  [1961]. 
These  authors  found  that  m  —  0.00016  for  the 
Gutenberg-Bullen  A  continental  model. 

The  curves  8099S  of  Kuo  et  al.  resulted  from 
the  application  of  the  Bolt  and  Dorman  spheric¬ 
ity  correction  to  the  Rayleigh  mode  phase  ve¬ 
locities  for  Dorman’s  oceanic  case  8099.  As  was 
previously  stated,  oceanic  phase  and  group  ve¬ 


locities  for  case  122  have  been  calculated  iiring 
a  program  for  Rayleigh  waves  on  a  spherical 
earth.  These  calculations  agree  with  the  resiilfs 
presented  by 'Kuo  et  al.  aA  curves:  £9998;  Thusr 
the  correction  for  the  effect  of  sphericity  oh 
Rayleigh  waves  is  nearly  identical  for  these:  ge¬ 
ologically  reasonable  tnbdels,  oceanic  and  Con¬ 
tinental. 

In  contrast  to  the  correcti&n  for  sphericity  for 
Rayleigh  waves,  which  is  the  saihe  for  many 
models,  the  corrections  for  sphericity  for  funda¬ 
mental  and  higher-mode  Love  waves  and  for  the 
shear  modes  of  the  Rayleigh  type  are  quite  de^ 
pendent  on  the  chCice  of  physical  parameters 
in  the  mantle.  Calculations  for  cases  122  for 
oceans,  CANSD  for  the  Canadian  shield,  and 
the  Gutenberg-Bullen  A  model  show  quite  dif¬ 
ferent  sphericity  corrections  for  the  fiinda- 
mental  Love  mode.  The  fundamental  Love  mode 
for  the  oceanic  model  develops  a  channel  wave, 
the  short-period  components  of  the  G  wave;  the 
fundamental  Love  mode  has  no  channel  wave 
for  these  continental  models. 

Future  ■problems.  Several  problems  remain 
to  be  solved  for  the  mantle  beneath  the  oceans. 
For  example,  the  calculated  phase  velocities  of 
oceanic  Rayleigh  waves  of  Kuo  et  cd.  [1962] 
are  a  few  hundredths  of  a  kilometer  per  second 
high  for  periods  greater  than  80  seconds.  Simi¬ 
larly,  the  calculated  phase  velocities  for  the  Gu¬ 
tenberg-Bullen  A  model  \Bolt  and  Dorman, 
1961]  are  slightly  above  the  data  for  periods 
greater  than  150  seconds. 

In  addition,  the  Calculated  Love  wave  ph^ 
velocities  in  Figure  1  of  this  paper  (spherical 
model)  are  also  several  hundredths  of  a  kilome¬ 
ter  per  second  too  high  for  periods  greater  than 
150  seconds.  As  a  consequence,  the  calculated 
group  velocities  are  somewhat  low  for  periods 
between  100  and  250  seconds.  Thus  both  Ray¬ 
leigh  and  Love  waves  indicate  that  small  changes 
in  either  velocity  or  density  must  be  made  in 
the  theoretical  models.  Preliminary  investiga¬ 
tions  show  that  either  shear  velocities  or  densi¬ 
ties  must  be  altered  at  depths  in  excess  of  2Q0 
km. 

Between  10  and  20  seconds,  the  calculated 
Love  wave  group  velocities  in  Figure  2  lie  above 
much  of  the  data.  Adjustments  of  shear  veloc¬ 
ity  in  the  crust  and  sediments  may  remove  part 
of  this  discrepancy  .  Changes  may  also  be  needed 
in  the  shape  Of  the  low-velocity  channel.  These 
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adjustments  will  probably  somewhat  alter  the 
shape  of  the  particle  amplitude  curves  for  the 
fundamental  Love  mode. 

OLservationai  data  for  higher  modes  were 
not  considered  in  the  present  Gt.udy.  These  modes 
should  be  observed  and  the  observations  ex¬ 
plained  by  the  model  used  to  fit  the  fundamental 
modes.  Finally,  we  note  that  S  wave  data,  still 
*  unavailable  for  oceaiuc  areas  at  epicentral  dis¬ 
tances  between  1“  aiid  17®,  and  between  17° 
and  34°,  would  be  of  great  help  in  future  studies. 
H  Summary  and  conclusions.  1 .  Love  wave  dis¬ 
persion  for  a  spherical  earth  has  been  computed 
for  a  model  of  the  suboceanic  crust  and  mantle, 
case  122.  This  model,  which  is  similar  to  case 
8099  of  Dorman  et  al.  [1960],  has  a  lovr-veloc- 
ity  channel. 

2.  Theoretical  Love  wave  velocities  for  case 
122  explain  phase  and  group  velocity  data  for 
oceanic  Love  waves  over  a  brood  range  of  pe¬ 
riods  to  within  a  few  hundredths  of  a  kilometer 
per  second. 

3.  An  oceanic  Rayleigh  wave  program  for  a 
spiw.i.al  eaith,  wricton  by  L.  '.’iop,  has  been 
used  to  confirm  the  sphencauv  voiiCCucn  phase 
and  group  velocities  for  Dorman’s  case  8099 
(presented  by  Kuo  et  al.  [i960]  as  case  8099S). 
These  theoretical  results  explain  the  oceanic 
Rayleigh  wave  data  of  Kuo  et  al.  to  within  a 
few  hundredths  of  a  kilometer  per  second.  The 
Love  and  Rayleigh  wave  results  can  be  further 
improved  by  minor  modifications  of  the  model. 

4.  The  effect  of  sphericity  on  Love  waves  is 
found  to  be  of  great  importance  and  to  be  de¬ 
pendent  on  the  model  under  consideration.  For 
the  oceans,  the  differences  between  the  phase 

f  velocities  which  result  from  calculations  for  flat¬ 
layered  and  spherical  earth  models  oxtcnc!  to 
periods  as  short  as  10  seconds. 

5.  Thus,  since  a  single  isotropic  distribution 
accounts  for  oceanic  Love  and  Rayleigh  wavs 
data,  the  discrepancy  between  calculations  for 
flat  layers  is  resolved  by  sphericity  without  the 
need  for  a  difference  in  SH  and  SF  velocities 
in  the  mantle  beneath  the  ceean.s. 

0.  Brune  and  Dorman  [1982]  found,  for  the 
mantle  beneath  the  Canadian  shield,  (1)  that  a 
low-velocity  zone  is  rcqtiired  and  (2)  that  no 
(iim.rjnce  between  SII  and  SF  velocities  is 
needed.  Th\ic  their  study  and  this  paper  sliow 
that  these  two  conclusions  are  ^p^dic-able  to  large 
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7.  For  periods  shorter  than  GO  seconds,  the  G 

wai'e,  which  is  the  fundamental  oceanic  liove 
mode  in  this  range  of  periods,  develops  a  large 
•sf-condary  maximum  of  particle  amplitude  in  tlie 
low-velocity  channel  and  may  be  considered 
a  channel  wave.  For  periods  between  60  and 
300  seconds,  the  free  surface  and  the  steep 
gradient  in  ehear  velocity  down  to  approximately 
900  km  form  the  wave  guide  of  the  G  wave. 
At  longer  periods  the  partieic  'uotion  penetvates 
the  high-speed  lower  mardir-,  phesft  and 

group  velocities  to  rise  sharply. 

8.  Calculated  higher-mode  Love  waves  and 
the  shear  modes  of  Rayleigh  type  also  exhibit 
secondary  maximums  in  the  lov.^-^'filf;city  chan¬ 
nel  when  their  phase  velocities  are  slightly 
greater  than  the  channel  velocity  and  their  ver¬ 
tical  wavelengths  are  ie,<w  than  the  ebannel  tbiclr- 
ness.  These  secondarj-  maximums  Kirill  be  impor¬ 
tant  in  the  consideration  of  problems  of  excita¬ 
tion  and  source  depth. 

9.  It  is  useful  to  compare  the  mantle  veloci¬ 
ties  determined  from  surface  wave  data  for 
oceanic  and  continental  areas.  Examples  include 
case  122,  case  £0G9  of  Dorman  et  al.  [1960],  and 
case  6EGHP1  of  Aki  and  Press  [1961]  for  the 
oceans  and,  for  the  contiiients,  the  Canadian 
shield  model,  CANSD,  of  Brune  and  Dorman 
[19G2]  and  the  Gutenberg-BuUen  A  model  dis¬ 
cussed  in  studies  of  the  Chilean  earthquake. 
When  this  comparison  is  made,  it  is  confirmed 
that  channel  velocities  are  smaller  beneath  the 
oceans.  Further  study  is  needed  to  determine 
whether  these  channel  velocities  usually  occur  at 
-shallower  depths  beneath  ocean  basins  than  be¬ 
neath  continent.s. 
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